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1 Introduction 

We consider in this paper the Cauchy problem for the evolution equation 

d t u + u(-A) a/2 u + V • (uVK * u) = 0, (1.1) 

u{x,0)=u o {x), (1.2) 

which describes swarming, collective motion and aggregation phenomena in 
biology and mechanics of continuous media. Here x G M. d , t > 0, and u = 
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u(x,t) > is either the population density of a species or the density of 
particles in a granular media. 

When v = 0, equation ( 11. ip can be considered as either a conservation 
law with a nonlocal (quadratic) nonlinearity or a transport equation with 
nonlocal velocity, and its character depends strongly on the properties of the 
kernel K. A classical choice for K is K(x) = e - ^' or, more generally, K is 
a radially symmetric function of r = \x\. Nonin creasing kernels correspond 
to the attraction of particles while nondecreasing ones are repulsive. Local 
and global existence of solutions to the inviscid equation ( 11.11) [y = 0) has 
been thoroughly studied in [13] under some additional hypotheses on the 
kernel, see also [HE]- In particular, kernels that are smooth (not singular) 
at the origin x = lead to the global in time existence of solutions, see e.g. 
[31 [13] . Mildly singular kernels (e.g. C 1 off the origin, like K(x) = e _ ' x ') may 
lead to blowup of solutions either in finite or infinite time [U El El UHl HH EES] • 
Strongly singular kernels like potential type (arising in chemotaxis theory, 
cf. [HIT]) 

K(x)=c\xf- d , (1.3) 

with 1 < (3 < d (so, in particular, the Newtonian potential kernel K(x) = 
Cd\x\ 2 ~ d , d > 3), usually lead to finite time blowup of "large" solutions, see 

[ZlEltlH]. 

Equations (ll.ip with fractional diffusion term [y > 0) have been intro- 
duced in the physical literature and studied in, e.g., [HI El El ELL], starting in 
nineties of the 20th century. The linear term in ( 11.11) is described by a frac- 
tional power of the Laplacian operator in M. d (or, more generally, by a Levy 
diffusion operator) defined in the Fourier variables by 

H(-^ a/2 m) = \z\ a m, (i-4) 

with < a < 2. 

When v > and if is a radially symmetric and nonincreasing function 
of r = \x\ with a mild singularity at r = 0, equation (11.11) then features 
a diffusive term which spreads the distribution of particles and a nonlinear 
drift term which concentrates it, thus acting in the opposite direction. The 



2 



fundamental question for the Cauchy problem fll.ll) - fll.2l) is to decide whether 
u governed by the competition between the nonlinear transport term and 
the linear dissipative term, can describe aggregation phenomena or not. Of 
course, the answer may depend on the (regularity and) size of initial data. 

Typical approaches to prove a finite time aggregation include an extension 
of the method of characteristics [21 [15] , the energy method (e.g. [H El [HI [15] ) , 
and the moment (or virial) method. The latter has been first applied to mean 
field models for self-gravitating particles and chemotaxis systems, [5], and 
recently in [HIT]. At this point, we mention that the characteristics method 
cannot obviously be applied in the presence of diffusion. 

Our aim in this paper is to present a simple virial type argument showing 
finite time blowup of a large class of solutions of f ll.ip - fll.2p if v > and 
< a < 1. The results we obtained are similar to those in [15] but we believe 
that our proofs are more direct and simpler. In addition, our assumptions 
on the initial data (see (12. 3p and (I2.4p below) and the kernel (see (12 .ip and 
(12.21) below) are less restrictive than those in [H] (radial symmetry and high 
localization on the initial data, the kernel K being a nonincreasing function of 
r = \x\) and [T3] (existence of exponential moments or even compact support 
assumption and K(x) = e - '^). In particular, the kernel K is allowed to have 
a repulsive part. 

The case of the strong dissipation 1 < a < 2 and nonlinearities with 
potential kernels (11.31) has been considered in [71 El [EE] where threshold con- 
ditions on the values of a, j3, al have been determined so that solutions can be 
either continued indefinitely in time, or they can blow up in a finite time for 
suitable initial data. But, for weakly singular kernels as the ones considered 
in this paper, the strong dissipation 1 < a < 2 prevents finite time blowup 
and global solutions exist, see [TT] and [HI Theorem 3]. 

Notation. 

The integrals with no integration limits are meant as J Rd .... Various posi- 
tive constants are denoted by C; sometimes the dependence of C on param- 
eters is written explicitly, e.g. C = C £ , C = C(e). 
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2 Main results 



We do not consider here local in time existence of solutions, their positivity 
and mass conservation properties since these topics have been discussed in 
detail in, e.g., [JJ El El EE]. One should note that the existence theory in 
[T4"l [To] is developed in the spirit of arguments used for conservation laws in 
|13j . i.e. without taking into account regularization effects of the diffusion 
while p], El El [UJ employed those effects in a significant way. 

Our results employ a crucial property of the gradient VK(x) of the con- 
volution kernel in (11 .11) . namely the fact that —x/\x\ is its homogeneous 
part near the origin. More precisely, we will use two sets of assumptions on 
the kernel K. There is a locally Lipschitz continuous function k such that 
K(x) = k(\x\) for x e R d , and 

(HI) either: there is K > such that 

-K <k'<0 and k r ■:= - sup k' > (2.1) 

(o,R) 

for each R > 0; 

(H2) or: k'(r) = —k[(r) + k' 2 (r) for r > and there are K > 0, K x > 0, 
K 2 > 0, and 5 G [0, 1) such that 

\k\r)\ < K (1 + r 6 ) , K x < k[{r) , k' 2 {r) < K 2 r s (2.2) 

for all r > 0. 

We will consider solutions which are even in x which is implied by the 
assumption that the initial condition u is even: 

it (x) = u (-x), x e M. d , (2.3) 

together with the radial symmetry of the kernel K and the uniqueness of 
solutions to (ll.ip - (ll.2l) . Moreover, we need that 

Mi := / |a;| u (x) dx < oo. (2.4) 
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As we have already remarked, the total mass is conserved during the evolution 

of Qnp-flOD 

/ u(x,t)dx = M := / u (x)dx. (2.5) 

Now we are in a position to state main results of this work. 

Theorem 2.1. Assume that v = 0. Consider a nonnegative and integrable 
initial condition u ^ satisfying Ii2.3\) and \2.1$ . 

(a) If k fulfills ^2. 1\) or ^2. 2\) with = 0, then the solution u to the Cauchy 

problem U.l\) - I[l~2]) ceases to exist in a finite time. 

(b) // k fulfills Ii2.2\) and Mi is sufficiently small, then the solution u to the 

Cauchy problem U.l\ )- I[l~2)) ceases to exist in a finite time. 

Let us emphasize here that the assumptions ( 12.11) and ( 12.21) with Ki = 
apply to two different classes of kernels k: indeed, k' is required to be bounded 
in the former but can vanish at infinity (in the sense that kr might decay 
to zero as R — > 00). The growth condition is less restrictive for the latter 
but k' is not allowed to vanish at infinity. We also point out here that we 
do not know whether the smallness of Mi is a necessary condition for finite 
time blowup to occur when k fulfills ( 12. 2p . 

It follows from Theorem 12.11 that, in the absence of diffusion and if the 
kernel is attractive (k' < 0), finite time blowup takes place for any nonzero 
initial data while a partially attractive kernel seems to require the initial data 
to be sufficiently concentrated for this phenomenon to occur. Since diffusion 
is expected to act also as a repulsive term, localization of the initial data 
seems also to be needed for finite time blowup when v > 0, even if the kernel 
is attractive. Indeed, we have the following result. 

Theorem 2.2. Consider a nonnegative and integrable initial condition uq ^ 
satisfying ( tff.gj) and \2.1$ . Assume that v > 0, < a < 1, and that k 
fulfills either Ii2.1\) or Ii2.2\) . If M is sufficiently large and Mi is sufficiently 
small, then the solution u to the Cauchy problem U.l\) - I[l~2}) ceases to exist 
in a finite time. 
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Observe that, besides localization of the initial data as in the partially 
repulsive case, Theorem 12.21 also requires the total mass M to be sufficiently 
large. This is due to the fact that, in the proof, it does not seem to be 
possible to balance the contribution from the diffusion with that from the 
drift term. 

In contrast to [Hj and pJ3 Theorem 12], our conditions on uq guaran- 
teeing finite time blowup do not require the L 1 -norm of u to be smaller (in 
a suitable sense) than that of u (K * u ). We also improve [T^l Theorem 8] 
where u Q is assumed to be in L 1 (M' i ; e 2 ' x ' dx). 



3 Virial inequalities 

For 7 £ (0, 1] and x £ M. d we define W 7 (x) = w 7 (\x\) with 

iw 7 (r) = i((l+r) 7 - 1), r>0. 

Evidently, W 1 > is a Lipschitz continuous function which will be used 
as a weight function. We list below some properties of and w 7 we will 
repeatedly use in the remainder of the paper. 

Lemma 3.1. Consider 7 e (0,1) and a £ (7,1). Then (-A) a / 2 W 7 £ 
L°°(M. d ). 

Proof. Recall that, for x £ M d , the Levy-Khintchine representation formula 
reads 

{ -^W l(x) = -C (d ,a) I iy 

jR d \y\ 

with 



27r (d+2a)/2r\( 2 - a )/2)' 
see, e.g., [TOj Theorem 1] or [T2j. Given x £ M d , we set g = g(\x\) 
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max {1, \x\}, and use the monotonicity and subadditivity of r i— > r 7 to obtain 

i -, a ^ t(x) i < 1/ (i + N + M);-(i + W d!/ 



C(d,a) |V ; tv /i - 7 



< - / 7(1 + M) 7 " 1 M^dy 



7 JB(o, e ) 

-I 

7 JR d \B(0,rf 

< C(d) (1 + (xl) 7 " 1 ^ + C(d) 



+- / v ' ' ' - ' ■ _'_TJ-. v ' ' ' u d v 



(1 + 


x|) 7 + | 


y 


7 - (1 + 


x|) 7 




y 





1 — a 7(a — 7) 



< C{d,a,i) Q^ a 



and the right-hand side of the above inequality is bounded since a > 7. □ 

Additional properties of w 7 are summarized in the next lemma. 

Lemma 3.2. Consider 7 £ (0,1]. For each e > i/iere exists a constant 
C e > sitc/i i/iai £/ie inequalities 

(1 — n> 7 (r)) < (1 — 7) u> 7 (r) and u> 7 (r) < — r 7 < £ + C £ w 7 (r) (3.1) 

/ioW /or all r >0. For 5 £ [0, 7) and R> 1 we have 

r s <^ /or r>R. (3.2) 

Proof. The first inequality in (13.11) follows from the observation that the 
function f(r) = (1 — 7) ((1 + r) 7 — 1) — 7 + 7(1 + r) 7_1 satisfies f'{r) = 
7(1 - 7 )(1 + r) 7 ~ 2 r > and /(0) = 0. 

The second inequality in (13.11) is clear for small r > and suitably large 
C = C e , as well as for large r ^> 1. 

Finally, if R > 1 and r > R, we have 



l+r 



! (l + r) 7 r 7 "" 5 



u, 7 (r) = / s 7 " 1 ds > r (1 + r) 7 " 1 > 1 ; > - a 



r 



2 ■ ' 

from which (13. 2p readily follows since 7 > 5. □ 

Next we derive an identity involving the moment J 7 of a nonnegative 
solution u of (11.11) defined by 



J 7 (t) := / W 7 (x)it(x,t)da; (3.3) 
7 



whenever it is meaningful (e.g., if u e L 1 ((0,T); (1 + \x\) dx)). 
Lemma 3.3. For each t > 0, we have 

^(t) = -uV(t) + Axit) + A 2 (t), (3.4) 

where 

V(t) := [ u(x,t)[(-A) a/2 W^](x)dx, 

At(t) := [[(w' y (\x\)-l)k'(\x-y\)^--^^-u(x,t)u(y,t)dxdy, 
JJ Fl \x-y\ 

M{t) := / / k'(\x-y\)^-- — — —u{x,t)u{y,t)dxdy. 
J J Fl \x-y\ 

Proof. The evolution of J 7 is governed by (11.11) so that 

^(t) = -v f W 7 (x) [{-A) a/2 u] (x,t)dx 
dt J^d 

+ w' y (\x\)-^-ru(x,t)-(VK*u)(x,t)dx 
JR d Fl 

t) [(-A) a/2 W 7 ] (x) dx 
w' {\x\)k'(\x - y|) A • |^ — i) w(y, t) dx dy 

Fl F ~ i/l 
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-i/D(t) 



> 7 (M) - - • ^ — —u{x,t)u(y,t) dxdy 

' fi f — y\ 

f f xx — y 

+ JJ k'(\x-y\)-—- ■ u(x,t)u(y,t)dxdy, 

whence lpTl)l . □ 

The next step is to find suitable upper bounds for P, Al, and Ai- Such 
an estimate for T> follows at once from Lemma 13.11 and (12.51) , and reads 

v(t) < IK-aj^xL / <x,t)dx, 

V(t) < C{d,j,a)M provided a E (7, 1). (3.5) 
Now we turn to A\ and A2, and first consider the case where k satisfies 
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Lemma 3.4. Assume that k fulfills A2.1\) . Then, for any R > 1, 

< (1 — 7) M K I y (t) , (3.6) 

Proof. We infer from (12. II) . (12.51) . and the first inequality in (13.11) that 
AW < JJ{l-w' 7 {\x\))\k'{\x-y\)\u{x,t)u{y,t)dxdy 

< (l-j)K // u(x,t)u(y,t)dxdy, 



whence ( 13.61) . Symmetrizing the double integral A2, we obtain 
AO) = ^ JJk'(\x-y\) (-^ - |^|J • ||— u(y,t)dxdy 

Since k! is nonpositive by ( 12. II) and 

1< If! + 1^ , (x,y)GM d xM d , (3.8' 
F - 2/| 

we deduce from ( 12 . 1 j) and ( 12.51) that, for any i? > 1, 
A(t) < ^ JJk'(\x-y\) U-j^-j^) M(x,t) M (y,t)dxdy 



< 



- / / fc'(|z-y|) (1-1^7-77) u(x,t)u(y,t)dxdy 

2 J B(0,R) J B(0,R) V Fl 1 2/1/ 



- ~~rr / / ( 1_ "rT'"TT) u(x,t)u(y,t)dxdy 

2 ys(o,ji) •/ -B(o,ji) V fI I 2/1/ 

+K-2R / ( / (l - A • r-| ] «(!/, f) dy ) dx 

jR d \jR d \B(0,R) V Fl 1 2/1/ / 



< f M 2_ ^u{x,t) dx' 2 

+k 2 r 1(1 u{y,t)dy)u{x,t)dx 
Jm. d \Jm. d \B(o,R) u]^[K) 
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since the evenness of x i— > u(x,t) warrants that 

/ ^- u(x,t) dx = for t>0. (3.9) 

The proof of Lemma 13.41 is then complete. □ 

We now derive the counterpart of Lemma l3~4l when k satisfies the weaker 
assumption (12.21) . Though the proof roughly proceeds along the same steps 
as that of Lemma 13.41 it is more complicated because some terms involving 
k\ and k<i have to be handled separately. 

Lemma 3.5. Assume that k fulfills ( BQ| ). Then, for any R > 1, 7 £ (8, 1) 

and e £ (0, 1), there is a constant C' £ depending only on e, 7 and 5 such that 

Ai(t) < K [(1-7) (2M + M R s + J 7 (i)) +2 M R 5 ~~<] I 7 (i) , (3.10) 

Mt) < M (^j + /{2c: ) / ^ ) -(f " 2A ' 2£ ) m2 ' (3 - n) 

Proof. On the one hand, we infer from (12. 2p and (13. ip that, for R > 1, 
Ai(t) < JJ(l-w^{\x\)) \k'(\x-y\)\ u(x,t)u(y,t)dxdy 

< K ff(l-w'(\x\)) (l + \x-y\ 5 ) u(x,t)u(y,t)dxdy 



< (1-7)^0 // w 7 (\x\) (1 + M 5 ) u(x,t)u(y,t)dxdy 



+K / (1 - w'^lxl)) \x\ d u(x,t) / u(y,t)dy) dx 

JB(0,R) \JR d J 

+K [ (l-w^(\x\)) \x\ 5 u(x,t) ( f u(y,t)dy] dx . 

JR d \B(0,R.) \JR d J 
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We next use (12.51) . (13.11) . (13.21) . and the property w'^ < 1 to obtain 

A 1 (t) < (l-7)^o U(l + (l + \y\y) u(y,t)dy^ I y (t) 
+(l-j)MK / Wy(\x\) \x\ 5 u(x,t)dx 

JB(0,R) 

2M K [ .. .. . . , 

H — d^v^a / w i\\ x \) u(x,t)dx 



Rl -S 

< (l- 7 ) K (2 M + 7 7 7 (t)) I 7 (t) + (1 - 7 ) K R s M I y (t) 

< K [(1-7) (2M + M R 5 + J 7 (t)) +2 M R s ^'] J 7 (t) , 



hence ( IXTUD . 

On the other hand, after the symmetrization of the double integral in A2, 
it follows from the positivity of k[ in ( 12.21) that, for R > 1, 

^ = 2 1 jj k '(\ x -y\) (]ff "iff) «(^*)«(y,*)dxdy 

= 5 fjk'{\x-y\) M±M^l-^-^ n(x,t)n(y,t)dxdy 

-f- + 4 ( 1 - rr ' FT J u(x,t)u(y,t)dxdy 



2 Jb{o,r) Jb{o,r) \ x — y\ V fI 1 2/ 1 

+^ / (7 fc'dx-yD^Mfi-^-^,] 
2 Jb(o,r) \Jri\b(o,r) k — 2/1 V fI Is/I . 



2 jR d \B{o,R) \jR d \ x ~ y\ V M \y\ 



< -\ I [ Wx-yfiM+Mfx-*.^) U (x,t)u(y,t)dxdy 
2 J b(o,r) J b(o,r) k — 2/1 V Fl 

+^ |y'^(|x-y|)M±M(l-| f .^ «(*,*) U (y,t)d*dy. 
Recalling ( 12.21) . (12.51) . (13.81) and ( 13. 9ft . we proceed as in the proof of Lemma l3~!^l 
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to conclude that 



B(0,R) JB(0,R) 

f I 

B(0,R) JB(0,R) 



k 'i(\ x -y\) t| z jr ( x ~ o ' iff I u ( x ^) u (v^) dx ^y 



> Ki I / ( 1_ |f["]f|) M <FF)w(2/,£)dxdy 



> K 1 M 2 -2K 1 / (/ u(y,t)dy J «(x,t)dx 

</R d \</R d VB(0,R) / 



r2 „ M If i 



Using once more (12.21) . (12.51) and the obvious bound 

<M+M^_ * .»)=(* _J»).£z» < 2 , ( „ )e ^, 

F-y| V fI M/ vfI M/ F-y| 



we find 



5 //^-^^fO-R'lif) «0MMv,f)d*di/ 



2 

< A' 2 |x — y| 5 £) £) dx dy 

< K 2 JJ {\x\ 5 + \y\ 5 ) u(x,t)u(y,t)dxdy 

< 2 M K 2 J \x\ 5 u(x,t)dx 

< 2M K 2 I e M + - ( ^— Ut) I , 



7 V l e 

whence (EUTl) . □ 



4 Finite time blowup 

4.1 The inviscid case v = 

Now we are ready to prove the first blowup result for the inviscid model 
(ll.li - p.2j) with v = 0, and begin with the case of a nonincreasing kernel K. 
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Proof of Theorem \2.1\ (a). We argue by contradiction and assume the solu- 
tion u of fll.lR — (jl.2p to be well-defined for all times. Combining (13.41) . ( 13.6ft 
and ( 13.71) . we end up with 

k 2R M 2 



^(t) < A,, R (I 7 (t)) 



M 



K-2R 

w 7 (R) 



/ 7 (t) 



for all R > 1 and 7 £ (0,1). Since A 7i r is a nondecreasing function, we 
realize that we have I 7 (t) < -f 7 (0) + A 7i ^(/ 7 (0)) t for t > as soon as 
A 7j r(/ 7 (0)) < 0. Then, of course, 7 7 attains zero at some finite time to which 
is impossible for nonnegative regular solutions to (jl.ip - (11.2p . a contradiction 
with the global existence. 

We next observe that we can always find 7 £ (1/2, 1) and R > 1 such 
that A 7vR (/ 7 (0)) < or equivalently 

k 2R M 



I-7) K + 



K 2 R 

w 7 (R) 



1,(0) < 



Indeed, if 7 £ (1/2, 1), we have u> 7 (r) < r for r > and w 7 (r) > \frj1 for 
r > 1. Therefore, choosing i? > 1 such that Mi < (M \/~R)/8 and then 
7 £ (1/2, 1) such that (1 - 7) K Mi < (k 2R M)/4, we realize that 



1-7) K + 



K 2 R 

w 7 (R) 



A(0) < 



2k 2R 



Mi < 



« 2 * M 



With this choice of R and 7, we have A 7 ^(/ 7 (0)) < and the proof is 



complete. 



□ 



Proof of Theorem \2.1\ (b). We again argue by contradiction and assume the 
solution u of (ll.ll) - (ll.2p to be well-defined for all times. Combining (13.41) . 
( 13.1 Op and ( 13.111) . we end up with 



d^ 



(t)<A 7 , R JUt)) for t>0, 



(4.1) 



where 



M 



;i- 7 ) K (2 + R s ) 



2K n 



R^~ s w^R) 



K2C' 



1 — 7) K z — [ — - — 2 K 2 e ) M 
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for all R > 1 and e 6 (0,1). As before, the inequality (14.1 p contradicts the 
global existence of nonnegative regular solutions to (jl.lj) - fll.2p as soon as 
A 7 ,_r j£ (/ 7 (0)) < 0. Since A 7i # i£ is an increasing function and / 7 (0) < Mi, we 
have A 7 ^ e (/ 7 (0)) < A 7vRie (M 1 ). Observing that an appropriate choice of e 
(sufficiently small) and R (sufficiently large) warrants A 7 r £ (0) < 0, we thus 
have A 7j # e (Mi) < provided Mi is small enough. Hence, for such a choice of 
e and R, finite time blowup of the solution to (jl.lft - fll.2p occurs as claimed. 

Finally, if K 2 = 0, we may argue as at the end of the proof of Theo- 
rem 12.11 (a) to show that, given any nonzero initial condition uq, we may 
find R large enough and 7 close to one such that A 7j ij ]£ (J 7 (0)) < 0, which 
completes the proof. □ 

Clearly, the only term in ( 14. ip that prevents Theorem 12.11 (b) from being 
valid for an arbitrary nonzero initial condition w is the term K 2 C' £ J 7 which 
cannot be made arbitrarily small by an appropriate choice of 7, R, and e. 
This term reflects the deviation of k from being decreasing, and thus the 
partially repulsive behaviour of k. 

4.2 The dissipative case v > 

The second result applies to solutions with suitably large initial data in the 
dissipative case: 

Proof of Theorem {KM Assume first that k fulfills ( 12. lft . We argue by 
contradiction and assume the solution u of fll.ip - fll.2ft to be well-defined for 
all times t > 0. Combining (13.4)) . (I3.5p . (13.6ft and (13.7ft . we end up with 



^(t) < A 7> *(/ 7 (*)) := M 



7) Kq 



w 1 (R) 
k 2R M 



I 7 (t) 



2 

for all R > 1 and 7 e (0,a). As before, the above inequality contradicts 
the global existence of nonnegative regular solutions to fll.lft - fll.2ft as soon 
as A 7jR (/ 7 (0)) < 0, the latter being true if A JtR (Mx) < 0. Fix 7 G (0, a) 
and R > 1 and assume that M > 4u C(<i, 7,0) / k 2 r. Then, A 7> ^(0) < 
-Mv C(d,j,a) < so that A 7|R (Mi) < if Mi is sufficiently small. 
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If k fulfills (O, the proof is similar and relies on flUD, (133]) . fl3~TU]) and 



In contrast to the proof of Theorem 12.11 (a), we cannot play with the 
parameter 7 in the proof of Theorem 12.21 when k fulfills (12. ip . Indeed, 7 is 
limited by the constraint 7 < a, and cannot be chosen arbitrarily close to 
one. This explains the necessity to have sufficiently localized initial data in 
the sense that Mi is required to be small enough. 
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